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ABSTRACT 

Using long-duration observations in the He II 304A passband of the Solar and 
Heliospheric Observatory (SOHO) Extreme-ultraviolet Imaging Telescope (EIT) we 
investigate the spatial and temporal appearance of impulsive intensity fluctuations in 
the pixel light curves. These passband intensity fluctuations come from plasma emit- 
ting in the chromosphere, transition region and lowest portions of the corona. We see 
that they are spatially tied to the supergranular scale and that their rate of occurrence is 
tied to the unsigned imbalance of the magnetic field in which they are observed. The 
signature of the fluctuations (in space and time) is consistent with their creation by 
magnetoconvection forced reconnection that is driven by the flow field in the high-/? 
plasma. The signature of the intensity fluctuations around an active region suggest that 
the bulk of the mass and energy supplied into the active region complex observed in 
the hotter coronal plasma is supplied by this process, dynamically forcing the looped 
structure from beneath. 

Subject headings: Sun: chromosphere - Sun:magnetic fields - Sun:transition region - 
Sun:corona 



1. Introduction 



Recent results (i McIntosh et al.ll2006l . |2007|) have provided detailed observational support for 
the ( widely held) hypo thesis that the relentless action of magnetoconvection-driven reconnection 
(e.g. JPriest et al.ll2002|) supplies the bulk of the energy to the solar atmosphere through ejecta that 
are intrinsically tied to supergranular s patial scales. In d iscussing thei r results, derived from Solar 
and Heliospheric Observatory (SOHO, Fleck et aDll995 ) observations, Mcintosh et al. ( 2007 ) have 
proposed that the small-scale eruptive phenomena observed are linked to the ubiquitous appearance 
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of "sp icules" at the solar limb in eclipse or coronagraph observations (e.g.. ISecchill 1 8771 : [Roberts 



1945|) . The work presented in this paper was influenced by a movie presented on the SOH O internet 
sitellfrom the Extreme-ultraviolet Imaging Telescope (EIT; Delaboudiniere et al. 1995 ) and shows 
the Sun "percolating" with passband intensity brightenings in the chrom osphere and transitio n 
region that occur on spatial scales, commensurate with those discussed by (IMcIntosh et al.ll2007h . 
This Letter (and a subsequent Paper that i s in preparation - Mc intosh, Judge & Gurman 2007) 
sets out to test the hypothesis presented by (IMcIntosh et al.ll2007|) with the EIT data used to build 
this movie and several others like it while offering an interesting, alternative, but logical look at a 
possible source of energy and mass input of active regions by the relentless forcing of emerging 
(and existing) magnetic flux by the convective flow of the high-/3 plasma. 

In the following section, we will discuss the time-series observations from the 304A (He II) 
passband of EIT that form the basis of our study. In Sect. 3 we will investigate the spatial and 
statistical dependence of the He II 304A intensity light-curves to identify (and count) individual 
"events" and study the variance of the intensity in each spatial pixel over an appropriate period 
of time. In Sect. 4 we will demonstrate that the appearance of spatial patterns in the analysis 
of the 304A observations are specifically tied to the plasmas magnetic environment, the chro- 
mospheric supergranular network and thus to the relentless action of magnetoconvection: emerg- 
ing, shuffling and destroyi ng small-scale magnetic flux elements that form the "Magnetic Carpet" 
(e.g, ISchrijver et al.l Il997h . We deduce, from the near ubiquity of the process, that its role is 
not restricted to loading mass and energy into the quiet sun and coronal holes, but is a signifi- 
cant source of the small-scal e dynamics that are re sponsible for mass and energy (with help from 



the mechanism discussed by IJefferies et al 



Aschwanden. Winebarger & Tsiklauril (|2007|) . 



2006|) into active regions; in large part the thesis of 



2. Observations 

The analysis presented in this Letter is based on the analysis of a long-duration sequence 
of EIT 304A passband imaging, starting 2006/03/18 20:36UT and ending 2006/03/21 23:48UT. 
The subset of the EIT sequence analyzed here (2006/03/19 00: 1 1 - 2006/03/19 23:59) has a fixed 
cadence of twelve minutes and has no missing data blocks and covers the approximate evolutionary 
lifetime of a supergranular cell (~24 hours). We assume that the bulk of the EIT 304A passband 
consists primarily of emission from He II formed in the upper c hromosphere or low solar transition 



region at a temperature of ~50,000 K ( Mazzotta et al. 



19981) . However, aside from the strong 



emission from the He II 303. 78A line there is a contribution to the bandpass emission in an active 
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region from the Si XI 303.31 A line and the blended emission of Mn XIV and Fe XV at 304. 86A 
(all three contaminant lines are fo rmed in the hotter corona at ~1.5MK) in the ratio 10::1::0.06 
( Brosius. Davila & Thomas 1998|) . The contribution of He II emission relative to the emission 
lines of Fe X/Fe XI (i.e., the contribution from the FIT 171 A passband in second spectral order at 
about ~1MK) is 25 to 1 in an active region and 10 to 1 in the quiet sun (|Brosius. Davila & Thomas 



19981). 



The 304A sequence is augmented by coronal context images from the FIT 1 95 A passband 



and fu ll disk line-of-sight magnetograms (Bu) from the Michelson Doppler Imager (|Scherrer et al. 



19951) closest to the start of the sequence, 01:12UT and 01:36UT respectively. Figure [B shows the 
context images and subfield pointing for the 2006 March FIT sequence. In panels A through D of 
each figure we show the FIT 195 A coronal image, the first FIT 304A image in the sequence, the 



MDI image and a smoothed image ((i?||)^_9n: lMcIntosh et al.ll2006|) providing information 
about the net-magnetic flux balance on the supergranular scale (i.e., the magnetogram is convolved 
with circular kernel of diameter 20Mm). We also show the 1200"x 1200" field of view studied 
(black square in panels A and B, white in panels C and D). The thin black contour in panel D 
shows where (-B||)r=2o=0 G, the magnetic "neutral" line occurs. 



The FIT data sample presented (reduced using the standard package discussed in the FIT 
User's Guide - |\protect||ittp://umbra.nascom. nasa.gov/eit/) shows a selected subfield of the 
FIT field-of-view that has been tracked, extracted and co-aligned (to within half of a pixel) us- 
ing the Interactive Data L anguage solar data mapping software in the SolarSoft data analysis tree 
Freeland & Handvl d 19981) . 



3. Data Analysis 

The analysis of the FIT 304A data presented is very straightforward, is easily repeatable 
and has two parallel threads to explore the temporal and spatial behavior of the observations. In 
the first thread, we analyze the intensity light curves at each spatial pixel in the de-rotated, co- 
aligned and trend (a ten timestep box-car smoothed profile of the intensity) removed (I{x, y, t) = 
I{x,y,t)— < I{x,y,t) >io) datacube to count the number of significant brightenings over the 
observing period. Fach brightening counts as one instance of an "event" and is defined as a 10% 
change in the trend-removed ratio A/// (= [i{x,y,t — 1) — i{x,y,t)]/ i{x,y,t)). The second 
thread investigates the spatial dependence of the brightenings through the pixel-by-pixel intensity 
distribution, f{I(x±5x,y±5y, t)), where Sx and Sy are each two pixels. Using a slightly extended 
spatial range ensures a higher signal-to-noise level in the time series (five times the number of 
points) and ensures that the moments of the intensity distribution, or a curve fit, can be accurately 
estimated with a high degree of reliability. The accumulation of data explicitly assumes that the 
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plasma in the adjoining pixels is (in some way) physically coupled to that in the pixel being studied. 

In Fig. |2] we take an example pixel from the 2006 March sequence to demonstrate the quan- 
tities used in the analysis. Panel A shows the intensity timeseries data for the pixel over 20 hours 
(triangles) and the ten pixel boxcar trend (gray solid line). In panel B we show the variation of 
A/// over the time period as well as the 10% increase (dot-dashed) line we use to define oc- 
curence an "event", the total number of events in that pixel is then computed as half the number 
of times the A/// profile crosses the 10% increase line. In panel C we show the distribution 
of intensities for the pixel range, a Gaussian fit (thick gray line) to the distribution and the four 
moments (mean, standard deviation, skewness and kurtosis) of the distribution. We note that the 
use of skewness and kurtosis in this analysis is very much depend ent on having a large enough 
number of sampl es to get an accurate intensity distribution (see, e.g. jAbramowitz & Stegunll 19721 : 



Press et al.ll 19921) . While they are not commonly used they serve a purpose in the present analysis. 



We use a simple definition of skewness: it measures whether outliers in the intensity distribution 
are below (negative) or above (postive) the distribution mean, skewing it to the left or the right re- 
spectively. Similarly, the kurtosis indicates whether the distribution is more or less sharply peaked 
than a Gaussian (a kurtosis of zero is perfectly Gaussian) while a fiat value has a negative kurtosis 
and very sharply peaked distributions result in a positive kurtosis. 



4. Results 

Applying the diagnostics discussed above to the EIT 304A timeseries sample we can develop 
an understanding of the relentless nature of magnetoconvection driven energy release in the transi- 
tion region and lower solar corona through the transient brightenings that it produces. In Figure [3] 
we present the diagnostic maps for the EIT sequence. Panels A through F show the mean EIT 304A 
intensity, the normalized width of the intensity distribution (width / mean intensity), the intensity 
distribution skewness, distribution kurtosis, the number of "events" and the (-B||)r=2o respectively. 
In each panel we overplot the EIT 195A iso-intensity contour at 150DN (e.g.. Fig. [l]\) while, 
in panel F, we show the magnetic neutral line as a thin black contour (e.g., Fig.djD). The online 
version of the Journal includes a version of this figure which includes the watershed segmentation 

that are derived from the mean 304A intensity over 
the period studied. The 500"x 400" rectangle around the active region (dashed rectangle) is shown 
for reference. 

Panel B shows the normalized width of the intensity fluctuation in each sequence. We see that 
the net width of the intensity distributions are almost uniformly reduced in the coronal holes on- 
disk (inside the EIT contours) where the mean intensity is also much smaller. There is a prevalent 
pattern of "rings" in this panel that appear to be of a size commensurate with the supergranular 



supergranular boundaries (e.g., Lin et al. II2OO3I 
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network, visible in panel A. In addition, there are positions of very large normalized widths ("hot 
spots"; at positions such as [-250", -325"], [-300",50"]) which are co-spatial with another set of 
rings that are clearly visible in panels C (skewness) and D (kurtosis). The skewness and kurtosis 
rings underlie the positions of EUV Bright Points (BPs) in the 195 A image (Mcintosh 2007 in 
preparation). Filament channels present their own particular signature in the normalized width 
maps (see, e.g., [-600:-200", -400:0"]) as nearly continuous lines ("snakes") of low distribution 
width in the center of the filament outlined with piecewise continuous regions of large distribution 
width. Active and plage regions (e.g., [0:200", -100:100"]) are most easily identified in panel E 
when the number of brightening events detected increases almost 75% above the mean background 
value (~35). 

Overall, we see that there is a striking correspondence between the number of events counted 
(panel E) and the var iation of the (-B||)r=2o (panel D). From the discussion in Sects. 4 and 5 of 
Mcintosh et al.l (|2007h we would expect this to be the case if magnetoconvection-driven reconnec- 
tion were responsible for their appearance; flux emerging into a region of larger imbalance would 
increase the probability of immediate reconnection. The same connection to patterning induced by 
magnetoconvection-driven reconnection can be attributed to the quiet sun distribution width rings. 
They occur almost exclusively at the interior of the supergranular boundaries and bear an uncanny 
resemblance to the loci of hig hest reconnection prob ability in the cartoon representation of the 
reconnection process (Fig 12 of iMcIntosh et al.ll2007[) . In this Letter we focus on the intensity di- 
agnostic patterns observed in the vicinity of an active region and offer a new insight into the mass 
and energy loading of active regions from their base, leaving other features for a future publication 
(Mcintosh, Judge & Gurman 2007 - in preparation). 

In what follows, we assume that the emission in the 304A bandpass, outside of the brightest 
active regions pixels, is dominated by H e II; we believe this assurn ption is justified both by the 
response of the FIT multilayer coatings (IDelaboudiniere et al.lll995h and spectroscopic measure- 
ments of quiet and active Sun intensities iBrosius. Davila & Thomasl(ll998b . Similarly, we assume 
that the constant driving of the emerging magnetic flux by the convective flow in the high plasma- /9 
will not switch off unless the flow itself is suppressed, i.e., in a sunspot spanning at least one or 
two supergranular cells in diameter. From additional context images supplied on the SOHO website 
(\protecthttp://sohow ww.nascom.nasa.gov/cgi-bin/get_sohoJmages?summary -1-200603 19 ) or on 
the Active Region Monitor 

(l\p r o t e c t tittp ://www. solarmonitor.org/index.php ?date=200603 19), there is a complex of very 
small sunspots underlying the active region, but that these spots are not of sufficient size to impede 
the flow significantly outside the spots themselves. 

Figure m shows the maps of the intensity distribution diagnostics in and around a relatively 
simple looking active region complex [NOAA ARs 10860 and 10862]. The panels of Fig. |4]show 
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the mean EIT 304A intensity, 195 A image, event location and the (-B||)r=2o maps for the active 
region. We see that spatial distribution of events (panel C) shows two lobes of intense activity 
where the number of events increases considerably (~60 events per pixel) over that of the quiet 
sun (~45 events per pixel) and the supergranular cell interiors in particular (~35 events per pixel). 
The event pattern is visibly correlated to the structure and amount of imbalance in the opposing 
polarity lobes of the active region itself with a measurable drop overlying the (-B||) r=2o neutral line 
(panel D) despite the bright emission over the same location in panels A and B. The location of 
the events tell us a great deal about how mass and energy are transported into the active region 
and plage complexes. This correspondence may also shed light on why neighboring coronal loops 
(anchored in the same active region) have apparently different mass loads and temperatures (e.g.. 



Klimchukll2006l and most TRACE observations). Panels E and F display enhanced distribution 



skewness and kurtosis in regions that bisect the major magnetic neutral line [100", 0"]. We specu- 
late that this, highly impulsive, but small (not a large number of events) magnitude intensity change 
in the bright emission above the sunspots has a significant contribution from the coronal emission 
contaminating the the EIT 304A passband (see the movies accompanying Fig.|3]), but acknowledge 
that information is really limited and this r nust be verified by model ing a nd more detailed obse rva- 
tions in the passband, such as discussed by Judge & Pietarilal i2004|); fpietarila & Judgel too4) and 
those from the EUNIS sounding rocket ( Thomas & DavilaboOll) . 



5. Discussion & Conclusion 



If we assume that we are indeed observing the impulsive mass loading and thermal heating 
caused by magnetoconvection-driven ejecta that are trapped belo w the coronal loop arca de (an ex- 



treme case of the quiet Sun heating and mass loading presented iMcIntosh et al. 



20071), then the 



most activity takes place in the largest regions of (-B||)r=2o- We would expect this to be the case 
when the emerging flux is almost immediately forced to reconnect in those locations because the 
net reconnection rate must increase as 1)^=20 increases. Similarly, the magnitude of the recon- 
nection event (the energy input into the system) will be greater at those locations, as there as there 
is a larger reservoir of magnetic imbalance to exploit, the drop in the event production rate around 
the neutral line to the quiet sun value appears to add weight to this argument. The increased vari- 
ation of the field inclination in the evolving high-/? m agnetic field will als o result in an increased 
"leakage" of p-modes in and around the active region (IJefferies et al.ll2006b that contain an energy 
flux may well dwarf that of the forced reconnection that initiated their propagation along the field 
lines. 



The net result of the footprint dominated energy release discussed here is that we will ob- 
serve coronal loops within the active region (and plage) complex where there is contrast in the 
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thermal and mass load based on the contrast in the magnetic flux available where they are an- 
chored. Further, we believe that the hotter (\ow-0 ) corona may well be sustained by an entirely 
different physical energy release (e.g., nanoflaring IParkeii 1 1988|) where the energy (and mass) is 
subsequently (and efficiently) redistributed by thermal conduction along the lines of force. In this 
case we have a "two stage" heating process for the active S un (as in the scenarios advocated by 
KlimchuklllOOd : lAschwanden. Wineba rger & Tsiklaurill2007b . There is a dominant hydrodynamic 
forcing of the plasma that initiates the mass and energy loading from the chromosphere and a 
secondary process that acts to redistribute the small percentage of the mass and energy along the 
magnetic fibers that compose the low-/5 portion of the active region. The second stage, intrinsi- 
cally tied to individual lines of fo rce and their interaction, should help resolve the apparent spat ial 
thickness of low-/5 coronal loops (|L6pez Fuentes. Klimchuk & Demoulinll2006l : lDeForestil2006() . 



Further, the results presented here may offer some perspective on the ano malously hig h 
UV/EUV helium emission observed on the Sun compared to theoretical models (|Jordanlll975|) . 
We note that the combination of magentoconvection-forced reconnecti on ejecta and leaking p - 
modes in a class of chromospheric spicule provide evidence in support of lJudge & Pietarilal (120041) 
that spicule forced neutral helium diffusion across the chromosp heric magnetic topology is the 
significant contributor to quiet Sun helium emission. In addition, ljudge & Pietarilal (|200 4') argue 
that the influence of the "photoionization-recombination" mechanism (IZirinlll975|) of overlying 
coronal radiation is negligible for the He II 304A emission that is formed in t he hyb rid physical 
region at the boundary of the chromosphere and transition region. iMauas et al.l (|2005h come to the 
same conclusion for active region UV/EUV helium emission. The other proposed helium emission 
enhancement mechan ism, "velocity redistr ibution" (from non-thermal plasma r notions in the line 
formation region, e.g. jAndretta et al. I boooh . is the result of "transient ionization" iPietarila & Judgd 
(|2004|) . We ask if it possible that this rapid transient ionization is a result of impulsively heating of 
the plasma as would happen locally in the forced reconnection ejecta we study in this paper? This 
question will be addressed at greater length in Mcintosh, ludge & Gurman (2007). 

To conclude. We have offered an interesting insight into the dynamic loading of (active region) 
coronal loop arcades from the high-/3 region in which they are tethered. In fact, we observe that the 
relentless loading is driven by the same mechanism responsible for the quiet Sun and coronal holes, 
while these have a clear topological difference on the global scale ("open" versus "closed" flux) 
an active region highlights the degree of energy input into the system from its footprint; we have 
seen that it scales as the degree of imbalance in the local magnetic field grows. Unfortunately, 
we cannot use the data presented in this Letter to discuss the exact micro-physics of how the 
dynamically forced loop footpoints in the active region transport their mass and energy into the 
structures observed in the hotter (low-/?) coronal plasma to answer, unambiguously, the question 
pose d in the title of this Letter, but th e detailed observations to be made by Hinode and FUNIS 
(e.g. jBrosius. Rabin & Thomasll2007r) will go a long way to resolving this picture. 
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Fig. 1. — Context images for the 2006 March 18 EIT 304A sequence. 
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Fig. 2. — Definition of the basic timeseries diagnostics used in this analysis (see text for details). 
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Fig. 3. — Distribution distribution maps for the March 2006 sequence. From left to right and top 
and bottom rows we show the mean intensity, normalized intensity distribution width, skewness, 
kurtosis, event location and (-B||)r=2o maps. In each panel we show the FIT 195A 150DN intensity 
level while, in panels F and F, we show the magnetic neutral line ({B^^)r=2o=OG). The dashed 
rectangle isolates the active region shown in Fig. |4l The online edition of the Journal has a version 
of this figure that shows the watershed segmentation boundaries derived from the 304A intensity 
image of panel A and also provides a movie of the entire 304A timeseries analyzed. The 500"x 
400" region around the active region (dashed rectangle) studied in this paper is shown for reference. 
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Fig. 4. — An expanded view of the active region present in the March 2006 sequence (cf. Fig. [3]). 



